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ABSTRACT 

Evaporation of water ice above 100 K in the inner few 100 AU of low-mass embedded protostars 
(the so-called hot core) should produce quiescent water vapor abundances of ^10~ 4 relative to H 2 . 
Observational evidence so far points at abundances of only a few 10~ 6 . However, these values are based 
on spherical models, which are known from interferometric studies to be inaccurate on the relevant 
spatial scales. Are hot cores really that much drier than expected, or are the low abundances an 
artifact of the inaccurate physical models? We present deep velocity-resolved Herschel-HIFI spectra 
of the 3i2-3 03 lines of H^O and H 48 (1097 GHz, E u /k = 249 K) in the low-mass Class protostar 
NGC1333 IRAS2A. A spherical radiative transfer model with a power-law density profile is unable to 
reproduce both the HIFI data and existing interferometric data on the H^O 3i3-22o line (203 GHz, 
E u /k — 204 K). Instead, the HIFI spectra likely show optically thick emission from a hot core with 
a radius of about 100 AU. The mass of the hot core is estimated from the C 18 J = 9-8 and 10-9 
lines. We derive a lower limit to the hot water abundance of 2 x 10 -5 , consistent with the theoretical 
predictions of ^10 -4 . The revised HDO/H2O abundance ratio is 1 x 10~ 3 , an order of magnitude 
lower than previously estimated. 

Subject headings: stars: formation — stars: protostars — circumstellar matter — techniques: spec- 
troscopic — astrochemistry 



1. INTRODUCTION 

In the embedded stage of low-mass star forma- 
tion, the central source is surrounded by a collaps- 
ing envelope t hat spans a wide r ange in densities and 
temperatures (A ndre et al.l I2000D . The inner region, 
known as the hot core o r hot corino (T > 100 K ; 
Walmslev fe Schilkd [1991 Ivan Dishoeck fe Blakd U99l 
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Ceccarellill200i r forms a crucial step in the flow of mat- 



ter from the cold outer enve lope toward the circumstel- 
lar disk (|Visser et al.l [2009) . Thermal evaporation of 
water ice at about 100 K shoul d produce a hot core 
water abundance of abou t 10~ 4 (jCeccarelli et al.ll"l996l ; 
iRodgers fc Cha rnlcv 2003) . Such abu ndances are indeed 
obser ved for high-mass protostars ([van der Tak et all 
120061: IChavarria et alj|2010t IHerpin et al.ll2012D . though 
not universally (jEmprechtinger et al.l I2013D . Abun- 
dances reporte d for low-mass hot co r es have so far not 
exceeded 10~ 5 (ICeccarelli et alll2000l:)Maret et all l2002l: 
iKristensen et al.l 120101: iCoutens et al.l I2012D . Are low- 
mass hot cores that much "drier" than expected? 

Measuring the hot core water abundance in a pro- 
tostar is not a trivial task, regardless of whether it 
is a low- or a high-mass source. Shoc ked gas tends 
to outshine the quiescent inner envelope (Mclnick ct al. 
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Coutens et al.l [20121 lEmprechtinger et al.l 120131 ). so one 
has to isolate the envelope emission from velocity- 
resolved spectra. Hot core abundances derived from 
spectrally unresolved data, su ch as from the In- 
frared Space Obs ervatory (ISO; ICeccarelli et al.l 120001 : 
IMaret et al.H2002T ). should therefore be treated with cau- 
tion. 

Another difficulty is the poorly known source structure 
on the spatial scales of the hot core: even if the column 
of quiescent hot water can be measured, there is no reli- 
able column of hot H2 to compare against. Abundances 
of water and other molecules are typically computed by 
way of spherical envelope models with simple density pro- 
files, con strained from single-d i sh dust continuum obser- 
vations (ICeccarelli et all 120001: Ivan der Tak et al.l l2000t 
Schoier et al.l 120021: IMaret et " all 120021 : ICoutens et al l 
20121) . Based on such a model, we previously de- 
rived a hot core abundance of < 10 -5 for the Class 
protostar NG C1333 IRAS2A (IKristensen et allllHOl: 
ILiu et all 1201 11) . However, spherical models do not re- 
cover the density enhancements measured with contin- 
uum interferometry i n the inner few 100 AU of various 
low-mass protosta rs (jJprgensen et all 12001 120051 120071: 
iChiang et alll2008T ) . The low water abundances may thus 
merely be an artefact of the inaccuracy of the adopted 
source models on the spatial scales of the hot core. 

Using interferometric observations of the 3i3~22o 
line of H^O at 203 GHz (E jk = 204 K) 
iJprgensen fc van Dishoecld (|2010bD and iPersson et all 
(2012) measured hot water column densities for the three 
Class protostars NGC1333 IRAS2A, 4A, and 4B. They 
also estimated H2 column densities from continuum in- 
terferometry on the same spatial scales to derive water 
abundances between a few 10 -9 and a few 10~ 6 . How- 
ever, these abundances are relative to the total amount 
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of gas in the inner ~100 AU, not relative to just the gas 
above 100 K. If a significant fraction of the material on 
these scales resides in an embedded disk or pseudo-disk, 
the bulk would be at lower temperatures and most wate r 
would be frozen out (jVisser et al.H2009t lllee et al.H201lD . 

In order to get additional constraints on the hot water 
in NGC1333 IRAS2A, we performed a 5.1-hr integra- 
tion of the 3i2-3 03 lines of H^O, Hf O, and H| 8 at 
1097 GHz with the Heterodyne Instrum ent for the Far- 
Infrared (HIFI; de Gr aauw et al.l 120101) on the Herschel 
Space Observatory (Pilbratt et al. 2010)0 This observa- 
tion is part of the key program "Water in star-forming re - 
gions with Herschel" (WISH; Ivan Dishoeck et al.ll20l"TI ). 
which aims to study the physics and chemistry of wa- 
ter during star formation across a range of masses and 
evolutionary stages. HIFI has high enough spectral reso- 
lution to disentangle the quiescent and shocked gas, and 
the deep integration allows for a detection of H^O. 

This paper presents the observations (Sections[5]and[3]) 
along with a detailed analysis (Section 2]). We develop a 
scenario in which the observed emission originates in the 
same gas as the 203 GHz line, namely in a hot core with a 
radius of about 100 AU. Our definition of a "hot core" is 
fairly loose and includes all quiescent (i.e., non-shocked) 
material hotter than 100 K, regardless of whether it is 
part of the envelope or an embedded disk-like structure. 
A simple spherical model fails at reproducing both data 
sets simultaneously, so we perform a model-independent 
analysis to constrain the distribution of water within the 
hot core. Comparing the water column density against 
Herschel-HIFI C 18 J = 9-8 and 10-9 data yields a 
lower limit on the water abundance of 2 x 10~ 5 relative 
to H2. As concluded in Section[5l the hot core in IRAS2A 
does not appear to be as dry as previously estimated from 
pure spherical models. 

2. OBSERVATIONS AND DATA REDUCTION 

NGC1333 IRAS2A is a low- mass Class protostar 
(L hol = 35.7 L , T bo i = 50 K; iKristensen etaTl [20121) 
located at a distance of 235 pc in Pe rseus at coordinate s 
3 h 28 m 55!6 by +31°14'37'.'l ('J2000: IBirota et al.l 120081) . 
Based on large-scale CO 3-2 maps, it is oriented close to 
edge on (ISandell et al.lll99l . IRAS2A was observed on 
2011 March 12 (obsid 1342215968) with HIFI on Herschel 
(jPilbratt et al.l r2010t Ide Graauw eTall [20101 ) in double 
beam switch mode with a nod of 3'. The primary tar- 
gets were the 3i2~3o3 lines of ortho-H^O and H^O at 
1097.365 and 1095.627 GHz in band 4 (E u /k = 249 K). 
The full spectrum covers a range of 4 GHz in each side- 
band, including the 3i2-3o3 and ln-Ooo lines of H3, 7 
and the 10-9 line of C 18 (Table |T]). The data were 
recorded in H and V polarization using both the wide- 
band spectrometer (WBS; 1.1 MHz resolution) and the 
high- resolution spectrometer (HRS; 0.5 MHz). The HIFI 
beam size at 1100 GHz is 19'/ 3 and the overa ll flux un- 
certainty is about 10% (jRoelfsema et al.|[2012D . 

The 3i2~3o3 line of H^O was previ ously detected in 
IRAS2A after a 30-min integration (|Kristensen et al] 
2010). The current observation had a total exposure time 
of 5.1 hr in order to detect the H^O line and to improve 

7 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with 
important participation from NASA. 
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Figure 1. Baseline-subtracted spectra of the 3i2-3o3 lines of 
H^O (bottom), H^ 8 (middle), and H^ 7 (top) in NGC1333 
IRAS2A. The H^O spectrum includes Gaussian profiles (green) 
fitted to the broad and medium components. Shown as an inset 
is a blow-up of the narrow feature in the H^O spectrum, with 
a Gaussian fitted in red. The narrow H^O feature is also fitted 
with a red Gaussian. The vertical dashed red lines mark the source 
velocity of 7.7 km s — 1 . The H^O spectrum is shown at the origi- 
nal spectral resolution of 0.14 km s . The other two spectra are 
rebinned to 1.0 km s — 1 and blown up by a factor of 10. 

the decomposition of the various velocity components in 
the H<!, 6 spectrum (Section [3]). The longer integration 
time also puts a more sensitive upper limit on the O 
line. 

The data were reduced with HIPE v8.2.1 (fOttl [20101 ) 
and exported to CLASS3 for analysis. The intensi- 
ties were converted from antenna to main-beam tem- 
perature scale through a main-beam efficiency of 0.74 
(jRoelfsema et al.l 120121 ) . The H and V spectra were av- 
eraged after individual inspection. The final data reduc- 
tion step was to subtract a linear baseline. 

3. RESULTS 

3.1. New Data 

The full HIFI spectrum as recorded with the WBS 
backend is shown in Figures [6] and [7] in the appendix. 
The strongest feature is the 3i2~3o3 line of H^O, peak- 
ing at 230 mK. The 3i2~3o3 spectrum for each water 
isotopolog is plotted in Figure [TJ Thanks to the excel- 
lent signal-to-noise ratio from the long integration (8 mK 

8 |http : //www . iram ■ f r /IRAMFR/GILDAS 
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Table 1 

Integrated Line Intensities (f T mh dv in K km s — 1 ) and 3<r Upper Limits for Narrow Emission in NGC1333 IRAS2A a 
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Telescope/ 
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E u /k 
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Observed 


Model 


Reference 13 






Instrument 
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| 


(eirCSGc) 


Intensity 


Intensity 




p-H-fO 


2()2-lll 


HIFI 


987.927 


100.8 


5.835 


-3) 


21.5 


< 0.309 


0.312 


1 


p-H-fO 


2n-2o2 


HIFI 


752.033 


136.9 


7.062 


-3) 


28.2 


< 0.583 


0.230 


1 


o-H^ 6 


3l2-221 


HIFI 


1153.127 


249.4 


2.634 


-3) 


18.4 


< 0.545 


0.352 


1 


o-H^O 


3l2~3o3 


HIFI 


1097.365 


249.4 


1.648 


-2) 


19.3 


0.342 ± 0.050 c 


0.342 


2 


p-H-fO 


331 -404 


HIFI 


1893.687 


410.4 


1.630 


-4) 


11.2 


< 0.531 


0.005 


3 


p-H^O 


2o2-lll 


HIFI 


994.675 


100.6 


6.020 


-3) 


21.3 


< 0.068 


0.005 


1 


tt1 8 i~\ 

P-H^ b O 


3l3-220 


PdBI 


203.408 


203.7 


4.812 


-6) 


0.83 


42 ± 8 


12 


4 


o-H^O 


3l2~3o3 


HIFI 


1095.627 


248.7 


1.621 


-2) 


19.4 


0.078 ± 0.011 


0.078 


2 


o-H^ 8 


4i4-3 2 i 


APEX 


390.608 


322.0 


3.143 


-5) 


16.0 


< 0.410 


0.205 d 


5 


o-H^ 8 


423-330 




489.054 


429.6 


6.887 


-5) 


0.83 




331 




o-HfO 


532-441 




692.079 


727.6 


1.478 


-4) 


0.83 




1.74 




P-H^ 7 


lii-Ooo 


HIFI 


1107.167 


53.1 


1.812 


-2) 


19.2 


< 0.037 


0.007 


2 


o-H^ 7 


3l2~3o3 


HIFI 


1096.414 


249.1 


1.635 


-2) 


19.3 


< 0.043 


0.005 


2 


c 18 o 


9-8 


HIFI 


987.560 


237.0 


7.380 


-5) 


21.5 


0.157 ±0.032 


0.043 6 


6 


c 18 o 


10-9 


HIFI 


1097.163 


289.7 


9.661 


-5) 


19.3 


0.150 ±0.050 


0.035 6 


2 



a Also listed are the frequency (u), upper-level energy (E u ), Einstein A coefficient (A u \), and beam diameter (8) for each line. 
The notat ion a (— b) denotes a X 10~ b . The model intensities are corrected for beam dilution and dust extinction as described 
in Section 14.31 

b 1: [Kristcriscn ct al. (20l3); 2: this work; 3: Benz et al. (in prep.); 4: [Persson ct al. (2012); 5: F. Wyrowski (priv. comm.); 
6: lYildiz et al.l [120101) 

c jT m t,di> for the broad and medium components is 3.34 ± 0.09 and 1.44 ± 0.08 K km s . 
d The intensity predicted for a fictitious 0'.'83 ALMA beam is 76.4 K km s~ 1 

e For an H2O/CO ratio of unity above 100 K. Most of the observed C 18 emission originates in gas below 100 K, which is 
not included in the model. 



rms in 0.5 km s _1 bins), the H^O spectrum now shows 
a distinct narrow emission feat ure on top of the medium 
and broad features detected bv iKristensen et"aT1 (j2010). 
The three components have full widths at half maximum 
(FWHMs) of 42.5 ± 0.8, 10.7 ± 0.4, and 3.6 ± 0.3 km 
s -1 . Table [T] lists the integrated intensities. The broad 
and medium FW HMs are the same as those reported by 
IKristensen et all for lower-excitation water lines. These 
components arise in shocked gas and are not considered 
here any further. 

The H^O spectrum shows a weak emission feature 
(7a) in both polarization filters, with the same intensity 
in the WBS and HRS r ecordings. The JPL and CDMS 
spect ral line databases (jPickett et ah 1998; Miiller et al.l 
[2001t ) show no other plausible features at this position in 
either the upper or the lower sideband. Figure Q] shows a 
Gaussian fit with an FWHM fixed at 3.6 km s _1 (taken 
from the H^O decomposition), because the feature is 
too weak for a reliable independent measurement of its 
width. The constrained fit has a peak intensity of 20 ± 4 
mK and an integrated intensity of 78 ± 11 mK km s _1 
(Table P. 

The H^O spectrum in Figure [1] shows an apparent 
emission feature blueshifted by 2 km s _1 from the source 
velocity, along with four features of similar strength at 
larger offsets. Three of them could represent lines of 
CH 2 DOH (E u /k = 229-294 K), which has b een detected 
befor e in single-dish spectra of IRAS2A (|Parise et al.l 
[20061: see also Figures M and [7]) . The JPL and CDMS 
databases offer no plausible identification for the fourth 
feature. We consider the 3i2~3o3 line undetected in 
H^O, but our analysis in Section |4] is unaffected by 
the question of weather it is real or not. The potential 



CH2DOH lines are not considered any further. 

The only other confirmed feature in the entire spec- 
tral setting is the C 18 10-9 line at 1097.163 GHz. It 
appears to be a blend of a narrow and a medium compo- 
nent, but the detection is too weak for an accurate de- 
composition. The intensity integrated over the full profile 
is 300 ± 30 mK km s _1 , split roughly halfway between 
the two components. Narrow and medium features have 
been detected p reviously in variou s low- and high- J CO 
isotopolog lines (jYildiz et al.H2010T ). The C 18 10-9 line 
is analyzed as part of a la r ger sa mple of lines and sources 
bv lSan Jose-Garcia et ail (|2013l ) and lYildiz eFaLl (|2013fl . 

The full spectrum in Figures |B] and [7] appears to contain 
a few more lines, but their nature is uncertain and not 
of interest to this paper. We merely mention two possi- 
ble identifications: H 2 CO at 1094.590 GHz (E u /k = 526 
K) and CH3OH at 1095.063 GHz (E u /k = 498 K). Both 
species have been detected previousl y in single-dish and 
interferometric studies of IRAS2A dMaret et al.l 120051: 
iParise et al.ll2006t Uorgensen et al.ll2007D . 

3.2. Complementary Data 

IRAS2A has b een targeted i n vario us other water lines 
with Herschel. iVisser et al.l (|2012f ) detected about a 
dozen spectrally unresolved lines with the PACS instru- 
ment (60-180 micron; E u /k up to 1750 K), attributed to 
shocked gas on 100-1000 AU scales and unrelated to the 
narrow emission analyzed here. No narrow emission has 
been detected in any of the water lines previously ob- 
served with HIFI IKristens en et al.l [20101 . Upper limits 
on such emission may be useful and are reported in Table 
[TJ Another HIFI upper limit is available for H^O 33i-4o4 
at 1894 GHz, observed alongside C+ 2 P 3/2 - 2 Pi/2 (Benz 
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et al. in prep.). In most cases, the upper limit in TableQ] 
is the 3(7 limit computed from the rms noise in 0.5 km s _1 
bins. This method does not work for the H 2 6 2o2-lii 
and 2n-2o2 lines, where the broad and medium emission 
features can hide a narrow feature stronger than three 
times the noise. Instead, we forced a three-Gaussian fit 
on these spectra, with the position and width of the nar- 
row feature fixed at the values found for the 3i2~3o3 line. 
The adopted upper limit is twice the integrated inten- 
sity of the force-fitted narrow component. In doing so, 
we assume the broad and medium components do not 
shield any narrow emission. The lio-loi, lii-Ooo, and 
2i2 _ loi lines show narrow a bsorption due to cold w ater 
vapor in the outer envelope (jKristensen et al.ll2010T > and 
are excluded from our analysis. 

Ground-based observations exist for two H 2 8 lines: 
3i 3 -2 20 at 203 GHz (E u /k = 204 K) and 4i 4 -3 2 i at 
391 GHz (E u /k = 322 K). The latter was undetected 
with APEX in a 16" beam, with a 3<r upper limit of 
410 mK km s -1 (Ta ble [H F. Wyrowski priv. comm.). 
IPersson et al.l (|2012D detected spatially resolved com- 
pact and extended emission in the 3i3-2 2 o line with the 
Plateau de Bure Interferometer (PdBI). The extended 
component is associated with the outflow and is ignored 
here. The compact component has an FWHM of 4.0±0.1 
km s~ x , comparable to the narrow component in the 
H 2 6 3i2~3o3 spectrum in Figure [TJ The compact PdBI 
component has an integrated intensity of 0.98 Jy km s _1 
or 42 K km s" 1 . Its spatial extent was marginally re- 
solved to a diameter of $.'83, prob ing material out to 
a radius of 100 AU. IPersson et al.l attributed the com- 
pact emission to a flattened inner envelope or pseudo-disk 
dominated by infall rather than rotation. 

4. ANALYSIS 
4.1. Spherical Model 

The HIFI 3i2~3o3 lines have upper-level energies of 249 
K. The bulk of the water in IRAS2A is present at tem- 
peratures below 100 K, where it generally exists as ice 
rather than vapor. Th e low gas-phase abundance in the 
cold envelope (~10 -8 ; IKristensen et al.ll2010t ILiu et all 
1201 lh . coupled with the low excitation temperatures, falls 
orders of magni tude shy of reproducing t he observed nar- 
row intensities (|van Kempen et a l. 2008). 

That leaves two options for the narrow 3i 2 -3q 3 emis - 
sion: a maser or the hot core. iFuruva et all ([2003) 
observed maser emission towards IRAS2A at 22 GHz, 
with intensities varying by an order of magnitude on 
timescales of about a year. The maser lines were seen 
sometimes at the source velocity and sometimes blue- 
shifted by a few km s~\ with FWHMs from 1.2 to 2.3 
km s -1 . The HIFI 3i 2 -3 03 and PdBI 3i 3 -2 20 spectra 
show no such blue-shifted emission and the emission at 
the source velocity is 2-3 times as broad as the 22 GHz 
maser lines. F urthermore, our two ep ochs of HIFI data 
(Figure [1] and IKristensen et al.l I2010D show no signs of 
variability on a one- year baseline. Hence, we conclude 
that the narrow 3i 2 -3q3 emissi o n is n ot associated with 
maser activity. IPersson et al.l ()2012D employed similar 
arguments to reach the same conclusion for their PdBI 
data. 

Instead, the narrow HIFI and PdBI lines are thermal 
and likely originate in the inner 100 AU of IRAS2A, in 



quiescent gas above 100 K. If the circumstellar mate- 
rial in IRAS2A is approximated as a spherical envelope 
with a po wer-law density profile, the 100 K radius lies 
at 94 AU (jKristensen et al.l 120121 ). While this matches 
the size of the emitting region seen with the PdBI, con- 
tinuum interferometry of IRAS2A has shown that the 
assumptions of spherical symmetry and a single power- 
law dens ity profile break down insid e a radius of l'/5 or 
350 AU (IJoreensen et al.ll200l 120051) . 

A simple test reveals that the spherical model also fails 
to reproduce the HIFI and PdBI data simultaneously. 
The test consists of synthesizing the 3i 2 -3o3 and 3i3~ 
220 spectra from a step abundance p ro file, as described 
in det ail by Ivan Kempen et al.l (|2008l ). IKristensen et al.l 
(|2012D constructed a source model for IRAS2A to fit the 
spectral energy distribution and sub-millimeter bright- 
ness profiles. The model envelope has a total mass of 5.1 
Mq, distributed along an r -1,7 power law. The temper- 
atures of the gas and dust are coupled and decrease from 
250 K at the inner edge (36 AU) to 10 K at the outer 
edge (17000 AU). The step abundance profile consists of 
a high abundance X m above 100 K (at 94 AU) and a 
lower abundance A out below 100 K. The HIFI and PdBI 
emission originate above 100 K, so the question we want 
to answer is whether the spherical envelope model can 
reproduce both data sets with a single Xj n . 

The molecular excitation and line emission are com- 
puted with the ID radiative tra nsfer code RATRAN 
(jHogerheiide fc van der Tak 2000). This is a standard 
approach for interpreting mol ecular line s pectra; see e.g. 
lYildiz et all (|2010l ). lLiu et al.l (|2011lh and lCoutens et al.l 
(l20ll for recent examples of C 18 0, HDO, and H 2 0. We 
set the gas-to-dust ratio to 100 and use OH5 opacities 
(|Qssenkopf fc Hennindl 19941 ). appropriate for dust grains 
with thin ice mantles. The ortho/para ratio of H 2 is 
thermalized and that of water is fixe d at 3 . Col lision 
rates are taken from IDubernet et ail (|2006l 120091 ) and 
IDaniel et al] (|2010l 1201 11) as compiled in the LAMDA 



database (jSchoier et al.ll2005D . The outer H^ 8 abun- 
dance does not affect the 3i 2 -3o3 and 3i3-22o lines and 
is fixed at 2 x 10 -11 , derived from abs orption seen with 
HIFI in lower rotational lines of water (jKristensen et al] 
[2011 ILiu etll][20rl . The inner H^O abundance is var- 
ied from 3 x 10" 10 to 3 x 10" 7 (2 x 10~ 7 to 2 x 10~ 4 for 
H 2 6 0) to cover the full range from a "dry" to a "wet" 
hot core. Finally, the synthetic spectra are convolved to 
the appropriate beam size and compared to the available 
observations. 

Figure [5] shows the convolved intensities for the H3; 8 
3i2 _ 3o3 and 3i3-2 2 o lines as function of the inner abun- 
dance. The horizontal dotted lines and bars mark the 
observed intensities and la uncertainties: red for 3i 2 - 
3o3 and black for 3i3~2 2 o. The abundances required to 
fit each line individually are marked by the vertical dot- 
ted lines and bars. Within the confines of the spherical 
envelope model, the H5, 8 3i 2 -3o3 spectrum from Figure 
Q] requires an inner H^O abundance ofl.2xl0~ 9 . How- 
ever, this underproduces the observed PdBI line inten- 
sity by a factor of 200. Matching the PdBI line requires 
a much higher inner H^O abundance of 2.1 x 10~ 7 , but 
then the 3i 2 -3o3 line comes out a factor of 6 too strong. 
The estimated uncertainty on the 3i 2 -3o3 intensity for 
H^ 8 is 17% (10% calibration, 14% statistical; Sections 
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Figure 2. Integrated intensities for the HIFI 3i2-3o3 line (red) 
and the PdBI 3i3-220 line (black) of H^O in the spherical envelope 
model of Kristcnscn ct al. (2012). The intensities, convolved to 
the appropriate beam size, are plotted as function of the hot core 
abundance of HJ^O (bottom axis) or H3j 6 (top axis) for a water 
o/p ratio of 3. The horizontal dotted lines and bars mark the 
observed intensities and lc uncertainties. The abundances required 
to fit each observation individually are indicated by the vertical 
dotted lines and bars. 

HandEU), so a factor of 6 is 35cr away. The RATRAN 
test clearly fails: a spherical envelope model with a sin- 
gle power-law density profile, such as used in previous 
attempts to derive hot core water abundances, cannot 
reproduce the single-dish and interferometric hot water 
observations simultaneously. 

4.2. LTE Analysis 

The mismatch of best-fit abundances to single-dish 
and interfero metric data has also be en seen for methanol 
in IRAS2A (|j0rgensen et all [2005h and for deuterated 
water in the Class protostar IRAS 16293 - 2422 
(jCoutens et al.ll2012t|Persson et al.ll2013l) . We thus have 
several pieces of evidence that simple spherical models 
break down on the spatial scales of the hot core. Until 
more appropriate source models can be generated - for 
example from sub-arcsecond PdBI and ALMA observa- 
tions of protostellar dust and gas - it is better to take 
a step back and see what we can learn about hot wa- 
ter in IRAS2A irrespective of the underlying density and 
temperature structure. 

4.2.1. PdBI Emission: Optically Thin 

In analyzing the PdBI data, iPersson et al.l (|2012l ) as- 
sumed optically thin emission in local thermodynamic 
equilibrium (LTE) and considered gas temperatures be- 
tween 50 and 250 K. The 3i3~22o line has a critical den- 
sity of a few 10 5 cm -3 . The critical density of the HIFI 
3i2~3o3 line is about a factor of 1000 higher, but with 
densities in the inn er envelope easily exceeding 10 9 cm~ 3 
(jVisser et al.ll2009T) , the assumption of LTE is probably 
justified and we adopt it here as well. 

The dust and gas temperatures are likely coupled at 
these high densities, so temperatures below 100 K can be 
excluded because of freeze-out. This also means that the 
0'.'83 emitting region resolved with the PdBI represents 
the entire reservoir of quiescent (i.e., non-shocked) hot 
gas in IRAS2A. The PdBI observations have a field of 
view of 25", a little larger than the HIFI beam of 19" at 
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1097 GHz. Hence, if any other hot gas had been present 
within the HIFI beam but outside the central 0"83, it 
would have been detected with the PdBI. 

Following iGoldsmith fc Langerl (|1999D . the optical 
depth of a molecular line is 

r» = ^^4(e h » /kT -l), (1) 

with v the frequency of the transition, A u \ the Einstein 
A coefficient, and Aw the observed FWHM. In LTE, the 
upper-level column density N u is 

N - Ngu c - E » /kT (2) 

with N the total column density, g u the degeneracy, and 
Q the partition function. The PdBI line is optically thin 
with t v ~ 0.1. For a water ortho/para ratio of 3, the ob- 
served PdBI intensity requires A^H^O) ~ 3 x 10 16 cm~ 2 
reg ardless of the exact temperature between 100 and 250 
K ()Persson et al.|[2012h . The 3 i2~3o3 lines have Einstein 
A coefficients a few thousand times faster than the 3i3- 
220 line (Table Q} and are optically thick for all three 
isotopologs. For example, at 100 K, r„ = 9.7 for H^O, 
34 for HfO, and 19 000 for H^ 6 0. The optical depths 
decrease by about 50% if the temperature is increased to 
250 K. 

4.2.2. HIFI Emission: Optically Thick 

If the HIFI emission in all three isotopologs is indeed 
optically thick and characterized by the same tempera- 
ture, the 3i2~3o3 lines would have had the same inten- 
sities. This is clearly inconsistent with the data (Figure 
[T|). However, the size of the optically thick region (ap- 
proximately the t u = 1 surface) may be different for each 
isotopolog, leading to different beam-filling factors. As 
argued below, this can explain the factor of 4.4 difference 
in observed intensities for H5, 6 and H^O as well as the 
non-detection of H.j> 7 0. 

Figure [3] shows a sketch of the flattened inner envelope 
of IRAS2A with the different emitting areas for the opti- 
cally thick 3i2-3o3 lines. The H^O 3i3~22o line observed 
with the PdBI is optically thin and traces the entire hot 
core to its outer radius of 100 AU. The extreme optical 
depth of 19 000 for H^O 3i2~3o3 means that its t v = 1 
surface, measured from the outside, essentially lies right 
at the hot core's outer edge. H^O therefore emits at 
close to 100 K. This is somewhat lower than the exci- 
tation temperature of 170 K assumed by IPersson et al.l 
(|2012D . but the differ ence has no substantial effect on 
our analysis. The temperature structure on these small 
scales is unknown, so we assume for simplicity that H^O 
and H^ 7 also emit at 100 K. The radius of the H^ 8 
r„ = 1 surface then has to be a factor of VL4~ = 2.1 
smaller (r « 47 AU) than that of H^ 6 (Figure [3]). The 
size difference would be smaller if we take a higher tem- 
perature for H;j, 8 0, but the available data do not allow 
for a more precise estimate. 

4.2.3. Water Distribution 

For an H^O column density of 3 x 10 16 cm -2 , the opti- 
cal depth in the H^ 8 3i 2 -3 03 line is 34 at 100 K (Equa- 
tion [T]). If the water is distributed uniformly throughout 
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Figure 3. Illustration of the hot core (T > 100 K) in IRAS2A. 
The red and blue areas show the t u = 1 surfaces for the 3i2— 3o3 
lines of H^O and H^O, with a difference in area that explains the 
different intensities observed with Herschel-HIFI (Figure \Q . The 
green area shows the maximum size of the H^O r„ = 1 surface 
allowed by the 3<r upper limit from HIFI. The red area is the same 
as that of the spatially resolved H3j 8 3i3— 220 emission detected 
with the PdBI (Persson ct al. 2012). The hot core is embedded 
inside a colder envelope (T < 100 K; shaded gray), extending out 
to 10 4 AU (not drawn to scale). 

the hot core, the t v = 1 surface for H^O would lie at 
l/34th of the way in from the outer edge at 100 AU, i.e., 
at r « 97 AU. The surface has to be at r « 47 AU 
to explain the observed H^O/H^O intensity ratio from 
Figure [IJ so the water cannot be distributed uniformly 
with radius. 

A power-law distribution, n(H20) oc r~ p , offers a sim- 
ple solution. This does not necessarily imply an abun- 
dance gradient, because it can be the tota l gas density 
itself that increases toward smaller radii (jVisser et al.l 
2009). We need to know the inner radius of the water 
reservoir in order to compute the power-law exponent 
p, but the inner radius cannot be measured from the 
PdBI data. A rough estimate would be 0.1 AU, the typ- 
ical inner edge of a ccretion disks around T Tauri stars 
(|Akeson et al.ll2005D . The requirement that only l/34th 
(3%) of the H^ 8 column lie at r > 47 AU then yields 
p = 1.2. The slope steepens to 1.4 if the inner radius is 
placed at 1 AU instead. 

Figure |4] shows the resulting optical depth profiles (t„ 
as function of r) for H^O and H^O. The r„ = 1 surface 
for H3, 6 (not plotted) lies at 100 AU, consistent with 
Figure [3] H^O reaches t u = 1 at 47 AU, as needed to fit 
the 3i2-3o3 observations. The t v = 1 surface for the least 
abundant isotopolog, H^O, ends up at 9.4 AU, providing 
110 times more beam dilution than experienced by H^O. 
This places the observable H^O 3i2-3o3 intensity about 
a factor of 10 below the 3cr upper limit from Figure [1] 
and Table HJ 

4.2.4. Summary of LTE Analysis 

In summary, the 3i2-3o3 emission observed with HIFI 
and the 3i3~22o emission observed with the PdBI can 
be explained simultaneously through an LTE analysis. 
The PdBI emission is spatially resolved and marks the 
outer boundary of the hot core in IRAS2A, with a ra- 
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Figure 4. Optical depth of H^O and H^O 3i2-3o3 measured 
inward from the outer edge of the hot core (Figure Q5). Both curves 
are for power-law water density profiles at an assumed emitting 
temperature of 100 K. The dotted lines mark the t u = 1 radii. 

dius of about 100 AU. The 3i3-22o line is optically 
thin and implies column densities of 2 x 10 19 cm~ 2 for 
RfO, 3 x 10 16 cm~ 2 for H j^O, and 8 x 10 15 cm" 2 for 
FP^O (|Persson et al.l [20121 ). The 1097 GHz HIFI lines 
are optically thick in all three isotopologs. The observed 
line intensities are nonetheless different because each iso- 
topolog becomes optically thick at a different radius, so 
that Hj O experiences less beam dilution than H^ 8 and 
H^ 7 (Figures 1 and 0J. 

The precise three-dimensional distribution of water 
within the hot core cannot be ascertained from the HIFI 
and PdBI data. The simplest solution is a radial power- 
law dependence, rt(H20) oc r~ p , with an exponent of 
about 1.2. We emphasize, however, that this power law 
is nothing more than a convenient approximation to the 
real 3D source structure in the hot core. The spheri- 
cal radiative transfer model from Section 14.11 fails not 
because its density power law in the inner 100 AU is a 
little too steep, but because the assumption of spheri- 
cal profiles for the density, temperature, and abundance 
does not hold on scales smaller than a few 100 AU. The 
implications of this conclusion are discussed in Section 
1431 



4.3. Dust Shielding and Observable Line Intensities 

The scenario from the previous section matches the line 
intensity ratios for the HIFI 3i2~3o3 isotopolog lines, but 
what about the absolute intensities of those and other 
lines? For optically thick emission in LTE, the intrinsic 
line intensity is the intensity of a blackbody at the gas 
temperature: 



■thick 



2fc^ 2 



B^T) . 



(3) 



The optically th in equation is taken from 

IGoldsmith fc Langerl (fl999h : 



thii 



hc 3 N u A 



ui 



(4) 



with all symbols as defined as in Equation [TJ The observ- 
able intensity T m b is the intrinsic intensity Tr divided by 
the relevant beam-filling factor. 
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For a gas temperature of 100 K, the optically thick 
H^O 3i2~3o3 emission has an integrated intensity of 320 
K km s _1 before applying beam dilution. The beam di- 
lution factor is (19.3/0.83) 2 = 540, reducing the intensity 
to 590 mK km s" 1 . This is still a factor of 1.7 stronger 
than observed (Table [IJ , which we attribute to extinc- 
tion by dust. The implied optical depth of the dust is 
r d = In 1.7 = 0.55 at 1097 GHz. 

Single-dish and interferometric continuum data of 
IRAS2A provide a lower and upper bound to the ob- 
served dust optical depth. The singl e-dish data, pro- 
cessed through the spherical model of iKristensen et all 
(2012), give a pencil-beam H2 column density of 3.7 x 
10 23 cm" 2 and r d = 0.26 at 1097 GHz. The interferomet- 
ric data reveal an embedded pseu do-disk with a gas mas s 
of 0.056 M Q within ~100 AU (|j0rgensen et al'l [2009h . 
which translates to r d = 2.3. The actual optical depth 
encountered by molecular line emission from the hot core 
has to lie between these two limits, so a value of 0.55 is 
quite reasonable. 

The final test is to check whether our hot core scenario 
also matches the upper limits on narrow emission in other 
water lines. Again assuming a gas temperature of 100 K, 
we repeat the above exercise to estimate intrinsic line in- 
tensities. Most lines are optically thick, so the location 
of their t v = 1 surfaces is calculated for the same r" 12 
power-law water distribution as in Figure HI The observ- 
able line fluxes, corrected for bea m dilution and dust ex- 
tincti on (based on OH5 opacities; Os senkopf fc Hennind 
fl99l . arc listed in Table [TJ The predicted intensity of 
the H^O 2 2-ln line at 988 GHz is right at the upper 
limit of how much narrow emission can be hidden in the 
observed broad and medium components. All other pre- 
dictions fall below the upper limits. Our scenario of a 100 
AU hot core with a power-law distribution of water (Fig- 
ures [3] and 0]) is therefore consistent with the combined 
PdBI, APEX, and Herschel data sets. 

The high optical depths for the water isotopolog lines 
accessible with HIFI are not unique to IRAS2A. For ex- 
ample, spatially resolved data from the Atacama Large 
Millimeter/submillimeter Array (ALMA) show that the 
column density of hot water in IRAS 16 293 — 2422 is 
a fac tor of 30 higher than in IRAS2A (jPersson et al.l 
120131). The HIFI w ater lines observed in this source 
(jCoutens et al.l 120121 ) therefore likely suffer from the 
same problem. A combined analysis of the ALMA and 
HIFI data is recommended for the best overall under- 
standing. 

The conclusions from this paper can be tested with 
ALMA, in particular regarding the excitation and spa- 
tial distribution of hot water. The 532-441 line of H3> 8 
at 692 GHz (E u = 728 K) is accessible in band 9 and 
has been observed in the hot core of IRAS 16293 — 2422 
(jPersson et al.l 120 13| ). Within a few years, receivers in 
band 8 (385-500 GHz) will be able to access the A14- 
3 2 i and 4 2 3-3 3 o lines of H^ 8 at 391 and 489 GHz 
(E u /k = 322 and 430 K). All three lines are predicted 
to be optically thin in IRAS2A and strong enough to be 
observable (Table [J). ALMA's spatial resolution of O'.'l 
or better in bands 8 and 9 will enable an investigation 
of the hot water in IRAS2A at an unprecedented level of 
detail. 
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Figure 5. H2O/CO column density ratio in the hot core re- 
quired to match the observed C ls O 9-8 and 10-9 line intensities if 
only a certain fraction of those intensities originates in gas hotter 
than 100 K. The shaded area brackets the fractions computed by 
Yildiz ct al. (2010, 2013). The dotted line marks the conservative 
lower limit to the H2O/CO ratio discussed in the text. 

4.4. Water Abundance and Deuteratium Fractionation 
in the Hot Core 

iPersson et al.l (|2012[ ) derived an H^O abundance of 
4 x 10~ 6 in the inner 100 AU of IRAS2A by comparing 
their H.j> 8 column density from the P dBI to the dust 
mass measured in a similar beam by Jorgcnsc n et al.l 
(2009). However, this is not the same as the hot core 
abundance: the simple picture painted in Figure [3] be- 
lies a more complicated three-dimensional structure, in 
which much of the gas and dust may be below 100 K. 
We present here an alternative approach to derive the 
real hot core water abundance, i.e., relative to only the 
gas above 100 K. 

H 2 cannot be observed directly, so CO is used as a 
proxy for the hot core mass. The best available tracers 
of the hot CO column are the optically thin 9-8 and 10- 
9 lines of C 18 {E n = 237 and 290 K), observed with 
HIFI in IRAS2A. There are two caveats, however: the 
hot CO abundance itself is uncertain by a factor of a 
few and some fraction of the 9-8 and 10-9 line intensi- 
ties origina tes at temperature below 100 K (jYildiz et al.l 
120101 12013D . We address the latter issue first. 

The 9-8 and 10-9 lines have critical densities of about 
10 6 cm -3 , so the assumption of LTE is justified. Invoking 
again an excitation temperature of 100 K, Equation [4] 
is used to predict the C 18 column density required to 
match a certain line intensity. In the hypothetical case 
that all of the observed 9-8 and 10-9 emission arises in 
the hot core (about 150 mK km s _1 for each line; Table 
CO), we get iV(C 18 0) 2 x 10 17 cm" 2 (corrected for beam 
dilution and dust extinction) or an H2O/CO abundance 
ratio of ^0.2. Figure [5] shows how the derived H 2 0/CO 
ratio changes if only a certain fraction of the 9-8 and 
10-9 emission arises in the hot core. For example, the 
H 2 0/CO ratio is 0.5 for a fraction of 50%. 

Full radiative transfer calculations on the spherical 
envelope model of IRAS2A show that about 10-30% 
of the observe d C 18 Q 9-8 and 1 -9 em ission originates 
above 100 K (jYildiz et al.l [20T1 [2T5T1 . This range is 
marked by the shaded area in Figure [5] and implies an 
H2O/CO abundance between 0.8 and 3.6. However, as 
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noted in Section 14.11 the spherical model is reliable only 
for radii larger tha n about 350 AU, where the dust tem- 
perature is ^45 K (Kristen sen et alj|2012l) . This can be 
used for a firm limit: the real fraction of 9-8 and 10- 
9 emission arising above 100 K cannot be larger than 
the fraction arising above 45 K in the spherical radiative 
transfer models. That latter fraction i s ~65-80% for the 

9- 8 l ine and ~90% for the 10-9 line (jYildiz et al.|[20Tol 
2013). According to Figure the lower limit on the hot 
core H2O/CO abundance ratio is then 0.25. 

lYildiz et al.l (j2010ft derived a C 18 abundance of 1.5 x 
10" 7 (8 x 10" 5 for CO) for all gas above 25 K, a factor 
of 3 below the carbon and oxygen elemental abundances. 
This value is well constrained for the bulk of the gas 
between 25 and 100 K, but remains somewhat uncertain 
for the hot core above 100 K: the hot core represents only 
a small mass fraction and contri butes little emission eve n 
to the C 18 9-8 and 10-9 lines (|Yildiz et al.l(2lMl2?)Tl . 
so varying the hot core CO abundance by a factor of a few 
would not affect the overall line intensities significantly. 

Taking the CO abundance of 8 x 10~ 5 in IRAS2A at 
face value, the results from Figure [5] yield a lower limit to 
the hot core water abundance of 2 x 10 -5 . This is a con- 
servative estimate, as the actual fraction of C 18 9-8 and 

10- 9 emission arising above 100 K is likely to be smaller 
than the fraction of 90% derived from the 45 K cutoff 
above. The available PdBI and HIFI observations are en- 
tirely consistent with a water abundance of (1-2) x 10~ 4 
in the hot core, as expected for the evapora tion of water 
ice above 100 K (jRodgers fc Charnleyl 20031) and as typi- 
cally found in high-mass protostars ( van der Tak et al.l 
[200a iChavarn'a etaTI [20101 iHerpin et all 12012ft . The 
uncertainties in the H2O/CO ratio and the hot core 
CO abundance also allow for water abundances down to 
about 10 -5 . Hot water abundan ces of less than 10~ 5 , as 
derived previou sly for IRAS2A ([Kristensen et al.l 120101 : 
ILiu et all 120111) or reported for other low-mass sources 
(iCeccarelli et al.ll200d iMaret et al.l 120021: iCoutens et al.l 
12012ft . are possible only if the inner CO abundance is 
significantly overpredicted by spherical models. 

Lastly, we revisit the question of deuterium frac tion- 
ation in the hot core of IRAS2A. ILiu et al.l (|201lD pre- 
sented five rotational lines of HDO measured with single- 
dish facilities. The best tracer of hot HDO is the 3i2~22i 
line at 226 GHz with an upper-level energy of 168 K. 
Once more invoking LTE at 100 K, the observed intensity 
of 0.50 K km s _1 translates to an HDO co lumn density 
of 2 x 10 16 cm" 2 (Equation |4|). ILiu et all derived a hot 
HDO/H2O ratio of > 0.01 based on a preliminary water 
abundance of < 10 -5 , which is now known to be too low. 
A more reliable deuterium fraction can be computed di- 
rectly from the water column density of iPersson et al.l 
(2012), bypassing the uncertainties in the H 2 0/CO ra- 
tio and the CO abundance. With AT(Hf O) = 2 x 10 19 , 
the new HDO/H2O abundance ratio in the hot core of 
IRAS2A is 1 x 10~ 3 . This is the same leve l of deuterium 
fracti onation as in IRAS 16293 - 2422 (jPersson et al.1 
12013ft and a factor o f 2 higher than the upper l imit fo r 
NGC1333 IRAS4B (IJorgensen fc van Dishoeckl l2010al) . 
both of which are based on intcrfcromctric data alone. 



4.5. Implications 



Continuum interferometry of IRAS2A and other low- 
mass embedded protostars has revealed the presence 
of a compact density enhancement, typically on ^100 
AU scales, that is n ot accounted for by simple spheri- 
cal envelope models (jJorgensen et al.ll2004L 120051 120071 
IChiang et al]l2008h . For a given column density of hot 
water, correcting for the higher density would lead to 
a lower water abundance than what is found from the 
spherical models. Our conclusion in the previous sec- 
tion is exactly opposite: the new hot water abundance, 
based on the combined PdBI and HIFI data for H3, 8 
and C 18 0, is higher than that obtained previously from 
single-dish data only. 

The reason for this apparent discrepancy is the tem- 
perature structure in the inner 100 AU of IRAS2A. With 
an H3, 8 column density of 3 x 10 16 cm -2 and an H3, 6 
abundance of 2 x 10~ 5 or higher, the column density of 
hot H 2 is 8 x 10 23 cm -2 or lower. The total gas col- 
umn within 100 AU, based on continuum interferome- 
try, is 4 x 10 24 cm -2 (jjprgensen et al.l 12009ft. Fr o m thi s 
we reach the same conclusion as IPersson et all (|2012ft : 
only a limited fraction of all qui escent materia l inside of 
100 AU is hotter than 100 K. IPersson et al.l derived a 
hot fraction of 4% for an assumed water abundance of 
1 x 10 -4 . Given the uncertainty in the water abundance, 
the actual hot fraction can be anything from 2 to 20%. 

The corresponding physical picture is not that of a 
spherical hot core, because in such a geometry it would 
be difficult to have the bulk of the material at temper- 
atures below 100 K, while still keeping some hotter gas 
out to 100 AU. An embedded disk or pseudo-disk is more 
likely, with a few per cent of its mass heated to above the 
evaporation temperature of water. This scenario is qual- 
itatively cojisisten^ 

mode ls dBrinch et all 120081 : 1 Visser et al.1l2009t lllee et all 
I201H lHarsono et al.ll2013ft . Specifically, it suggests that 
the bulk of the disk starts cold, with little water ice evap- 
orati ng as material acc retes from the larger-scale enve- 
lope (|visseLeEaD[200l). 

Regarding the chemistry, the lower limit on the hot 
core water abundance is consistent with the value of 
~10 -4 expected from the ev aporation of water ice above 
100 K. iStauber et all (|2006l ) presented models of hot and 
cold envelope chemistry in the presence of UV and X-ray 
fields. For an X-ray flux of 10~ 3 erg s _1 cm -2 (equiva- 
lent to an isotropic X-ray luminosity of 3 x 10 28 erg s _1 
at 100 AU), the hot water abundance is reduced from 
10~ 4 to 10" 6 on a timescale of 10 4 yr. However, the gas 
above 100 K is essentially in freefall, leading to a dynam- 
ical timescale of p erhaps only 100 yr C Sc hoier et al.ll2002l : 
iVisser et al1 l20091 . That is still long enough for water to 
evaporate from the grains, but not for it to be destroyed 
by X-rays. Our relatively high hot water abundance can 
be sustained by continuous replenishment of icy grains 
from further out in the envelope. 

Several complex organic species have been observed 
in IRAS2A, includi n g CH3OH, CHcjCN , and CH3OCH3 
(IMaret et all 120051: iParise et all 120061: iBottinelli et all 
120071 : Uorgensen et all 120071: IPersson et al.l 12012ft . The 
line widths of a few km s -1 and the excitation tem- 
peratures of about 100 K mark IRAS2A as a proto- 
typical hot core or hot corino. This is in contrast to 
other protostars such as NGC1333 IRAS4A and 4B, 
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where the observed complex organics appear to be as- 
socia ted with outflow-induced shocks ([J0rgensen et al.l 
2007). The aforementioned dynamical timescale of 100 yr 
in the hot core of IRAS2A is too short for the formation 
of complex organics in the gas phase after evaporation of 
the icy grain mantles. Instead, the organics observed in 
IRAS2A are likely to be so-cal led first-generation species 
(|Herbst fc van Dishoeckl [20091 ) . formed via grain-surface 
chemistry in the cold outer envelope and released into 
the gas phase in the hot core. 



5. CONCLUSIONS 

This paper presents Herschel-HIFI spectra of the 3i2~ 
3 03 lines of H^ 6 0, H^ 7 0, and H^ 8 (1097 GHz; E u /k = 
249 K) in the low-mass Class protostar NGC1333 
IRAS2A. The deep 5.1-hr integration brings out a nar- 
row emission feature (FWHM = 3.6 km s _1 ) in the H3, 6 
spectrum on top of two broader components. The H3, 8 
spectrum shows a narrow feature only and the line is not 
detected in H3, 7 0. The broad components seen in H3, 6 
are associated with shocked gas and are not discussed 
here further. 

Guided by spatially res olved H^O 3i3-2 2n emission 
detected with the PdBI (jPersson et all 12012( 1. we at- 
tribute the narrow emission from the HIFI H^O and 
H^O spectra to a hot core (T > 100 K) with a ra- 
dius of about 100 AU. A spherical envelope model with 
a power-law density profile and a step abundance for 
water is unable to reproduce the HIFI and PdBI in- 
tensities simultaneously. This is consistent with previ- 
ous c omparisons between single-d ish and interferometric 
data ([Jprgensen et al.ll2004l . l2005t l and is additional proof 
that these simple models cannot be used for an accurate 
analysis of protostellar hot cores. 

We present an alternative scenario where the HIFI lines 
are optically thick and the difference in observed intensi- 
ties is due to different areas of the r„ = 1 surfaces of the 
three isotopologs. The hot core is moderately shielded 
by dust with an optical depth of 0.55 at 1097 GHz. The 
beam-diluted and extinction-corrected emission from the 
hot core matches the intensities of the PdBI line, the two 
detected 3i2-3o3 HIFI lines, as well as upper limits on 
various other lines targeted with HIFI. 

Because the detected HIFI emission is optically thick, 
the column density of water in the hot core (100 AU 
radius) is unchanged from the values of 3 x 10 16 cm" 2 for 
ffifO or 2 x 10 19 cm -2 for H^O derived bv lPersson et al.1 
(2012). We derive CO column densities above 100 K from 



C 18 J = 9-8 and 10-9, and find a conservative lower 
limit to the H2O/CO abundance ratio of 0.25. The hot 
CO abundance is 8 x 10 -5 (jYildiz et al.l 12010ft . so the 
water abundance in the hot core can be anywhere from 
2 x 10~ 5 to 2 x 10" 4 . Abundances of less than 10~ 5 , 
as previously derived from a spherical model, are only 
possible if the CO abundance above 100 K is significantly 
lower than that between 25 and 100 K. The revised water 
abundance implies a hot HDO/H2O ratio of 1 x 10~ 3 , 
an order of magnitude lower than earlier estimates and 
similar to interferometric HDO/H2O measurements in 
other Class protostars. 
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APPENDIX 
FULL SPECTRUM 

The full spectrum from our Herschel-HIFI observations of NGC1333 IRAS2A covers frequencies from 1094 to 1098 
GHz in the lower sideband (LSB) and 1106 to 1110 GHz in the upper sideband (USB). With a total integration time of 
5.1 hr, this is one of the deepest spectra obtained with Herschel for a low-mass embedded protostar. The full spectrum, 
rebinned to a velocity resolution of 1.0 km s , is shown in Figures [6] and [71 The molecular line frequencies from Table 
[T]are red-shifted by the source velocity of 7.7 km s _1 and marked with dotted lines. Also marked are the positions of 
six other transitions located within the two sidebands: one of H2CO, one of CH3OH, and four of CH2DOH. 
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Figure 6. Full baseline-subtracted HIFI spectrum for IRAS2A at 1094-1098 GHz (LSB, bottom axis) and 1106-1110 GHz (USB, upper 
axis) at 1.0 km s — 1 resolution. The dotted lines mark the locations of the transitions from Table[T]and of some organic species, all corrected 
for a source velocity of 7.7 km s — 1 . The spectrum itself is plotted at the original observed frequencies. 
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Figure 7. Same as Figure [6] but with the spectrum broken up into four chunks of 1 GHz to zoom in on the weaker features. 
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